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Abstract. If all stars form in clusters and both the stars and the clusters follow a power law 
distribution which favours the creation of low mass objects, then the numerous low mass clusters 
will be deficient in high mass stars. Therefore, the mass function of stars, integrated over the 
whole galaxy (the Integrated Galactic Initial Mass Function, IGIMF) will be steeper at the high 
mass end than the underlying IMF of the stars. We show how the steepness of the IGIMF depends 
on the sampling method and on the assumptions made for the star cluster mass function. We 
also investigate the 0-star content, integrated photometry and chemical enrichment of galaxies 
that result from several IGIMFs, as compared to more standard IMFs. 
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1. Introduction 



In recent years, a series of papers (Kroupa fc Weidner 2003 Weidner fc Kroupa 2004 



Weidner fc Kroupa 2005 and Weidner fc Kroupa 2006 WK06 from now on) have pro- 



posed that the stellar content of an entire galaxy may not be well described by the same 
initial mass function (IMF) that describes the distribution of stellar masses in the star 
clusters, where these stars form. The reason is that star clusters also form with a cluster 
mass function (CMF), which is a power law with a power law index of —2. If the lower 
cluster mass limit is very low (in the mass range that is also occupied by single stars, 
i.e. below 100 Mq, and under the assumption that the most massive star in a cluster 
cannot be more massive than its host cluster), the low mass clusters will be deficient in 
high mass stars. Therefore, if the stellar content of all clusters is added up, resulting in 
the Integrated Galactic Initial Mass Function (IGIMF), the distribution of stellar masses 
may be steeper at the high mass end, depending on the exact shape of the CMF. 

The exact form of the IGIMF could have profound implications for the integrated prop- 
erties of galaxies. High mass stars are important for galaxies, both in terms of their energy 
output (massive stars are bright and explode in supernova explosions) and in terms of the 
chemical enrichment. Assuming an IGIMF instead of an IMF, the she supernova rate and 
chemical enrichment of galaxies are studied by Goodwin fc Pagel (2005)[ the relation be- 



tween Ha and UV luminosity in star forming disks by Pflamm-Altenburg fc Kroupa (2008) 



Recchi et al. (2009)) investigates [a/ Fe] abundance ratios for IGIMFs. 

In most of these earlier studies, the assumed CMF and sampling method are kept the 
same. For sampling they use the method 'sorted sampling', which was put forward by 
WK06, as they found it to best fit the relation between the most massive star in a cluster 
and the cluster mass. The CMF was assumed to be a pure power-law with index -2.2 and 
a lower cluster mass limit of 5Mq . 
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In this contribution, we will vary the sampling method and the CMF parameters and 
investigate the resulting IGIMFs. Also, we will estimate the observability of the effects 
of sampling method and CMF variations through the number count of 0-stars as will 
be observed by GAIA, the integrated photometry of distant galaxies and the chemical 
enrichment of galaxies (in terms of total gas phase metallicity). 

2. Constructing IGIMFs 

2.1. Sampling method 

As a default choice for the CMF, we will use the same CMF as used before: 

dN/dMci oc MJ^^'^, with a lower mass limit of 5Mq. For computational simplicity we 

use the Salpeter IMF for stars, between 0.1 and lOOM©. 

We test several sampling methods (italic parts are the descriptions used in the figure 
later on): 

(a) Random sampling until the cluster mass is reached. As this will never get exactly 
to the cluster mass aimed for, we cither include the last star {stop after), exclude it (stop 
before), only exclude it if that brings the total mass in stars closer to the cluster mass 
aimed for than including it {stop nearest) or excluding the last star (which first surpassed 
the cluster mass) at a 50% probability {stop 50/50). 

{b) Sampling a specific number of stars, estimated from the total cluster mass and 
average mass in the IMF (with the maximum stellar mass limited to the cluster mass or 
not {number and number unlimited)). 

(c) The ^sorted sampling^ method of WK06. 

{d) Analytic sampling, in which we assume that the IMF in a cluster is always sampled 
analytically, excluding stochastic effects. 



For more detailed descriptions, see Haas & Anders (2009) We use a Monte Carlo 



approach in order to make the numerical scatter smaller than the difference between 
IGIMFs. We investigate the dependence of the results on the underlying IMF (we test 
for Kroupa) and the results are virtually indistinguishable. 

2.2. Cluster mass function 

As the main cause of the lack of high mass stars are the low mass clusters, the behaviour 
of the CMF at the low mass end is very important. The default choice is an extrapolation 
of the observed CMF at higher masses {Md >^ lO'^Mo). down to the mass of the smallest 
observed star forming clumps in the Taurus- Auriga region. We therefore test, besides the 
default choice given above, also a lower limit of IM0 and SOMq for the clusters, and 
power law indices of -1.8, -3.2 and -4.2 (effectively changing the number of low mass 
clusters, compared to high mass clusters). 



3. Resulting IGIMFs 

In Fig. [1] we show the resulting IGIMFs from different sampling methods (left) and 
different CMFs (right). The value of the IGIMF is divided by the value for the underlying 
IMF, in order to enhance the visibility of the differences. Also, dotted lines indicate the 
slopes the different IGIMFs would have. 

It can be clearly seen that the different sampling methods give different IGIMFs. The 
recovered CMFs after filling the clusters with stars are indistinguishable from the input 
CMF. Depending on the method, the steep end has power law indices of about -2.6, 
as compared to the input -2.35. The 'sorted sampling' and 'analytic sampling' get even 
slightly steeper and deviate strongly from a power law shape. 

From the Monte Carlo realisations using different CMFs. it can be seen that the slope 
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Figure 1. The resulting IGIMFs (normalized to the input IMF) for different sampling methods 
(left) and CMFs (right). The indications correspond to the descriptions in the text. The dotted 
lines with indices indicate the steepness that line would have in an IGIMF. 

of the CMF determines the high mass end slope of the IGIMF, whereas the minimum 
cluster mass sets the stellar mass at which the steepening sets in. If clusters do not 
typically go down in mass all the way to the mass range of stars, sampling effects in 
the IGIMF will be negligible. The low mass end of the CMF is badly constrained. The 
IGIMF could be obtained observationally, giving an indirect measure of the low mass 
end of the CMF (under the assumption of a sampling method and an underlying IMF). 

4. Galaxy properties from IGIMFs 

Using the IGIMFs as IMFs in galaxy evolution models we are able to investigate the 
integrated properties of galaxies, if their total stellar initial mass distribution is given 
by an IGIMF rather than a more standard IMF. The high mass end of the stellar mass 
distribution is important for galaxies, as their light in many pass bands is dominated 
by the most massive stars for a wide range of star formation histories, their chemical 
composition depends on the metals expelled by the massive stars and the ISM is strongly 
influenced by the energy output of the massive stars exploding as supernovac. 

4.1. The number of 0-stars in the Galaxy 

As the most pronounced effect appears for the most massive stars, we first estimate the 
difference in massive star content of the galaxies. Specifically, we estimate the number of 
0-stars (with masses > ITM©) that GAIA will be able to observe. Assuming that GAIA 
will observe 10% of all 0-stars in the Galaxy, that they live for 10 Myr and that in the 
past 10 Myr the average SFR of the Galaxy was 1 MQ/yr, GAIA will be able to rule 
out several of the very extreme IGIMFs (from e.g. steep CMFs and sorted sampling) 
with very high significance. Judging between the several different sampling methods is 
difficult, as the differences are of the order of la, for purely poissonian errors on the 
number counts. 

If our assumptions regarding the observed fraction of 0-stars, the SFR and or the 
lifetimes of 0-stars are inaccurate, the resultant number of 0-stars varies linearly with 
either of them. The errors on the numbers scale with the square root of the numbers. As 
the numbers we calculate are of the order of 1000 to 2000, and the differences an order 
of magnitude smaller, the combined effect of our assumptions may be off by a factor of 
a few and the differences between all IGIMFs will still be of more than 2a significance. 

4.2. Integrated photometry of galaxies 

Using the GALEV evolutionary synthesis models (jBicker et al. 20041 IKotuUa et al. 2009| 
we are able to assess observable properties of galaxies, in which the stellar initial mass 
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distribution is given by IGIMFs, instead of more popular IMFs. The GALEV models 
follow the chemical enrichment and photometric evolution of galaxies self-consistently in 
closed-box models of galaxy evolution. They arc able to match a wide range of observed 
properties of galaxies of all Hubble types, like stellar mass, star formation and multi-band 
photometry. 

The ingredients of the models are star formation histories which depend on the gas 
mass present (no in- and outflows are assumed, nor are they necessary) and Hubble Type, 
spectral synthesis, and yield tables. Subsequent generations of stars are formed from the 
gas enriched by earlier generations. 

For all Hubble types, the difference in integrated photometry in any band is smaller 
than the observed galaxy-to-galaxy scatter as present in the HyperLeda databas^ (jPaturel et al. 2003)) 
Therefore, integrated photometry of galaxies will not be a discriminant between IMFs 
and/or IGIMFs. 

4.3. Chemical enrichment 

Many heavy elements arc mainly produced in supernova explosions of Type II (includ- 
ing Ib,c) and in the winds of massive AGB stars. With a deficiency in massive stars, 
as expected from the IGIMF, less production of metals is to be expected. The GALEV 
models do not follow the chemical enrichment on an clement by element basis, but rather 
follow the total mctallicity. Element abundances can then be obtained by assuming solar 
abundance ratios. 

With IGIMFs as determined above, metallicities of galaxies vary by half a dex at given 
gas mass fraction of the galaxies. For the IGIMFs resulting from extremely steep CMFs, 
the difference can be up to 2 orders of magnitude. At given gas mass fractions, the Hy- 
perLeda database gives a scatter of about half a dex as well. Therefore, the most extreme 
CMFs can be ruled out by the observed gas phase metallicities. In order to differentiate 
between several IGIMFs more precise measurements of metallicities are necessary. The 
difference between sampling methods is comparable to, or slightly larger than the spread 
that results from using individual metallicity measurements for the same galaxy, and is 
thus a possible discriminant between IGIMFs, provided that absolute determinations of 
gas phase metallicities are reliable. 
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